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New Results on the Hadronic Contributions 

TO a(M|) AND TO {g - 2)^ 



We reevaluate the dispersion integrals of the leading order hadronic contributions to the running 
of the QED fine structure constant a{s) at s = M|, and to the anomalous magnetic moments of 
the muon and the electron. Finite-energy QCD sum rule techniques complete the data from e~^e~ 
annihilation and r decays at low energy and at the cc threshold. Global quark-hadron duality is 
assumed in order to resolve the integrals using the Operator Product Expansion wherever it is 
applicable. We obtain Aahad(M|) = (276.3 ± 1.6) x 10"^ yielding a-i(M|) = 128.933 ± 0.021, 
and a}^^'^ = (692.4 it 6.2) x 10"^'^ with which we find for the complete Standard Model prediction 
= (11659159.6 it 6.7) x 10^^'^. For the electron, the hadronic contribution reads a^^'^ = 
(187.5 ± 1.8) X lO"^''. The following formulae express our results on the running of a at M| as 
a function of the input value for as{M^) and its error: 
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Introduction 



The running of the QED fine structure constant a{s) and the anomalous magnetic moment 
of the muon are observables whose theoretical precisions are limited by second order loop 
effects from hadronic vacuum polarization. Both magnitudes are related via dispersion 
relations to the hadronic production rate R in e~^e~ annihilation. While far from quark 
thresholds and at sufficiently high center-of-mass energy y/s, R{s) can be predicted by 
perturbative QCD, theory may fail when resonances occur, i. e., local quark-hadron duality 
is broken. Fortunately, one can circumvent this drawback by using e~^e~ annihilation data 
for R{s) and, as proposed in Ref. [|I|, hadronic r decays benefitting from the conserved 
vector current (CVC). With help from a moment analysis using weighted integrals over 
low-energy e~^e~ cross sections we showed in Ref. that the Operator Product Expansion 
(OPE) (also called SVZ approach 0]) can safely be applied down to energies of = 
1.8 GeV. It turned out that at this energy nonperturbative contributions to the dispersion 
integrals are very small. 

In this letter, we improve our previous determinations by using finite-energy QCD 
sum rule techniques in order to access theoretically energy regions where perturbative 
QCD fails. This idea was recently advocated in Ref. [Q. In principle, the method uses no 
additional assumptions beyond those applied in Ref. 0. However, parts of the dispersion 
integrals evaluated at low-energy and the cc threshold are obtained from values of the 
Adler D-function itself, where we assume local quark-hadron duality to hold. We therefore 
perform an evaluation at rather high energies (3 GeV for u, d, s quarks and 15 GeV 
for the c quark contribution) to suppress deviations from the local duality assumption 
by nonperturbative phenomena. We present a thorough evaluation of the associated 
theoretical and experimental uncertainties. 

Running of the QED Fine Structure Constant 

The running of the electromagnetic fine structure constant a{s) is governed by the renor- 
malized vacuum polarization function, 11^ (s). For the spin 1 photon, n^(s) is given 
by the Fourier transform of the time-ordered product of the electromagnetic currents 



j^^{s) in the vacuum (gV - ^V'') n^(g^) = i I d'^xe''i''{0\T{j^^{x)j^^{0))\0). With 



Aa{s) = —4:iiaRe[Il^{s) —11^(0)] and Aa{s) = Aaiep(s) + Aahad("S), which subdivides 
the running contributions into a leptonic and a hadronic part, one has 



where 47ra(0) is the square of the electron charge in the long- wavelength Thomson limit. 

For the case of interest, s = M|, the leptonic contribution at three-loop order has 
recently been calculated to be p 




a(0) 



(1) 



1 - Aaiop(s) - Aahad(s) 



Aaiep(M|) = 314.97686 x 10 
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Using analyticity and unitarity, the dispersion integral for the contribution from the 
hadronic vacuum polarization reads 



Aa.^(M|) = Re/ ds j^j^^,— . (3) 



and, employing the identity l/{x' — x — ie)e^o = P{l/(x' — x)} + in6{x' — x), the above 
integral is evaluated using the principle value integration technique. 



Muon Magnetic Anomaly 



It is convenient to separate the Standard Model prediction for the anomalous magnetic 
moment of the muon, = [g — 2)^/2, into its different contributions, 

«r = + + «r^' ' (4) 

where aj^™ = (11 658 470.6 ± 0.2) x 10^^° is the pure electromagnetic contribution (see fg] 
and references therein), aj^'^*^ is the contribution from hadronic vacuum polarization, and 
0'^'^^^ = (15.1 ± 0.4) X 10~^° P, H, |lOl accounts for corrections due to exchange of the 
weak interacting bosons up to two loops. 

Equivalently to Aahad(^z)) by virtue of the analyticity of the vacuum polarization 
correlator, the contribution of the hadronic vacuum polarization to can be calculated 



via the dispersion integral |TI 



ar = ^^^fiW. (5) 

4m2 



where K{s) denotes the QED kernel |T2 



K{s) = x^ (l-^] + {1 + x)2 + ln(l + x)-x + ^] + 7^^a;2 Inx , (6) 



2 ' ' \ x^J \ ' ' 2 1 - X 



with X = (1 — /3^)/(l + /3^) and (3^ = {1 — 4m^/s)^/^. The function K{s) decreases 
monotonically with increasing s. It gives a strong weight to the low energy part of 
the integral (|^). About 92% of the total contribution to is accumulated at cm. 
energies y/s below 1.8 GeV and 72% of a^"^ is covered by the two-pion final state which is 
dominated by the p(770) resonance. Data from vector hadronic r decays published by the 
ALEPH Collaboration |jl3| provide a precise spectrum for the two-pion final state as well 
as new input for the lesser known four-pion final states. This new information improves 
significantly the precision of the determination [|l|. 

Theoretical Prediction of R{s) 

Using a method based on weighted integrals over the low-energy e^e~ cross sections, we 
showed in Ref. that perturbative QCD is safely apphcable for the evaluation of the 
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dispersion integrals and at energies above 1.8 GeV, leaving out the cc threshold 
region where resonances occur. We stress that this approach uses global quark-hadron 
duality only. Deviations from local duality are averaged to a negligible contribution when 
calculating the integral. However, a systematic uncertainty is introduced through the 
cut at explicitly 1.8 GeV, i.e., non- vanishing oscillations could give rise to a bias after 
integration. In order to estimate the associated systematic error, we fitted different 
oscillating curves to the (rather imprecise) data around the cut region and obtained the 
following estimates 

A(Aahad(M|)) = 0.15xl0-\ 

Aal""^ = 0.24 X 10"^° , (7) 

from the comparison of the integral over the oscillating simulated data to perturbative 
QCD. These numbers are added as systematic uncertainties to the corresponding low- 
energy integrals. 

In asymptotic energy regions we use the formulae of Ref. |T^ which include mass 
corrections up to order to evaluate the perturbative prediction of R{s) entering into 
the integrals (H) and (|). 

Theoretical Approach Using Low-Energy Spectral Moments 

Following the suggestion made in Ref. we can write the following identity: 

F= j dsR{s)[f{s)-p.^^m{s)] + — j ^[PnA^o) - PnA^)]Duds{s) , (8) 
Ami \s\=so 

with P„,„(s) = !^dtpn,m{t) and f{s) = a{QfK{s)/{?,'K^s) for F = a^^^^^^, as well 
as f{s) = a(0)M|/(37rs(s - M|)) for F = A«had (M|) p^^, ^j. The contour integral 
runs counter-clockwise around the circle from s = sq — ie to s = sq + ie. The regular 
functions Pn,m{s) approximate the kernel f{s) in order to reduce the contribution of the 
first integral in Eq. (|D which has a singularity at s = and is thus evaluated using 
experimental data. The second integral in Eq. (H) can be evaluated theoretically apply- 
ing QCD. To suppress unphysical subtractions we use the Adler D-function, defined as 
D/.(s) = —12n'^sdllf.{s)/ds and RfXs) = 127rlmn/-(s + ie) for the set of quark flavours 
fi, to calculate theoretically the second integral in Eq. (|]), rather than the correlator 11/.. 



Using the OPE, the D-function for massive quarks is given by [|T^, [Tg, [T^l 




^(6 + 28^ + (295.1 - 12.3 n,) (^''' 
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+ 127r^i^ + IGTT^i^ \ , (9) 



S3 ' s4 ^ 

where additional logarithms occur when fi"^ ^ s and fi being the renormalization scaleQ. 
The coefficients of the perturbative part read do = 1, di = 1.9857 — 0.1153 rij, d2 = 
d2 + /5o^^/48 with /5o = ll — 2n//3 and nj the number of involved quark flavours, and 
d2 = -6.6368 - 1.2001 Uf - 0.0052 nj - 1.2395 (E/ Qff/Nc E/ Q/- The nonperturbative 
operators in Eq. @ are the gluon condensate, {{as/'n')GG), and the quark condensates, 
{mfqfqf). The latter obey approximately the PCAC relations 

{iTLu + nid) {uu + dd) ~ -2f^ml , 

ms{ss) ~ -/^(m|-m2), (10) 

with the pion decay constant /^r = (92.4 ib 0.26) MeV |T^. In the chiral limit the relations 
fn = fx and (uu) = (dd) = (ss) hold. The complete dimension D = 6 and D = 8 opera- 
tors are parametrized phenomenologically using the vacuum expectation values (Oq) and 
(Os), respectively. 

The functions Pn,m{s) introduced in Eq. (||) are chosen in order to reduce the uncer- 
tainty of the data integral. This approximately coincides with a low residual value of the 
integral, i.e., a good approximation of the integration kernel f{s) by the Pn,m{s). We use 

where the form (1 — s/sq) ensures a vanishing integrand at the crossing of the positive 
real axis where the validity of the OPE is questioned. Polynomials of order involve 
leading order nonperturbative contributions of dimension D = 2{n + 1). We therefore 
restrict the analysis to n < 2. Including additional powers n > 2 would require to fit 
the expectation values of the corresponding dimension D > Q operators as they are not 
known theoretically and a direct application of the results on r data ||2^ is dangerous 
since deviations from the vacuum saturation hypothesis used could be s-dependent. The 
new information is then spoiled since the additional fitted parameters deteriorate the 
accuracy of Aahad(^z) '^Ji^*^- Figure ^ shows the integration kernel of a^^'^ and two 
of its Pn,m approximations obtained from a minimization. At low energies pi,2o shows 
a clearly better agreement with the kernel than o so that the uncertainty from the 
data integral in Eq. (|]) is reduced. However, setting m > increases dramatically the 
theoretical uncertainties from as, ms and the nonperturbative operators which are then 
much larger than the errors of the pure data dispersion integrals. Therefore, we choose 
m = in the following analysis. 

Several polynomials are used as spectral moments in order to adjust a,|i'^pm^, 

^ The negative energy-squared in D(—s) of Eq. (^ is introduced when continuing the Adler function 
from the spacehke Euclidean space, where it is originally defined, to the timelike Minkowski space by 
virtue of its analyticity property. 
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Figure 1: The integration kernel K{s)/s of Eq. (^) and approximations according to the 
definitions Eqs. (pH]) and ([M])- 



Aahad(M|)[2„ 



as well as the gluon condensate simultaneously in a combined fit 



using Eq (P). The experimental and theoretical correlations between the moments are 
calculated analytically from the total experimental covariance matrices and the variation 
of the theoretical parameters within their uncertainties. 



Theoretical Approach Using Dispersion Relations 

Another closely related method is the approximation of the integrals (^) and (|^) via the 
dispersion relation of the Adler D-function 



RAs) 



D,iQ^) = Q'Jdsj^ 



Ami 



+ Q 



2\2 



:i2) 



for space-like = —q^ and the quark flavour /. Using the above integrand as approxi- 
mation of the integration kernels in Eqs. and (|), Eq. (H) becomes 



so 

F = JdsR 

Ami 



Data/ 



AfQ' 



s + Q 



2\2 



+ Af lDuds{Q^)-Q^Jds 



QCD/ 

uds \ 



so 



[s + Q 



2\2 



, (13) 



with a normalization constant Ap to be optimized for both F = a^^fzm^, ^ = 

Aahad(^z)[2m^, y^]- The last integral in Eq ( |13|) corrects the contribution from the D- 
function and is calculated theoretically. One notices that in contrast to Eq. (j^), where 
only global quark-hadron duality is assumed, Eq. ([I3| ) requires local duality to hold at 
Q^. One therefore needs to choose large against the scale where nonperturbative 
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phenomena appear. On the other hand, too large deteriorates the approximation of 
the kernel f{s) in Eq. (|1^) enhancing the experimental uncertainty on F. 

It is interesting to test whether derivatives of the Adler function (|12[), defined as 

^'"<«'' ^ '"'-"'^Qy^'' = . (14) 

Am) 

with Tn{s,Q'^) = {Q'^Y^^ /{s + Q'^Y^'^, ameliorate the evaluation of the hadronic contri- 
butions Higher derivatives > 15 indeed improve the approximation of the kernel 
function f{s) by ApTn{Q'^) as can be seen from Fig. |l], and thus reduce the experimental 
uncertainty from the first integral on the r.h.s. of Eq. (|T3p. Unfortunately, it increases 
the theoretical errors of the D-function in Eq. (|13[). In order to demonstrate the effect 
we use the theoretical prediction (||) which, setting q;s(M|) = 0.1200 and a/Q^ = 3 GeV, 
yields for m, d, s quark flavours the moments 

Mo(9 GeV^) = 2.190 + Q.2Sml + 0.15((a,/7r)G'G) + 1.8(m,gg) + 0.3(06) + 0.05(08) , 
M4(9 GeV^) = 0.468 + 0.45m, + 0.46((a,/7r)GG) + h^m^qq) + 2.3(06) + 0.7(08) • 

Taking as uncertainties: Aa,(M|) = 0.002, Am, = 0.07 GeV, A((a,/7r)GG) = 0.02 GeV^ 
A{mgqq) = 0.3 X 10"^ GeV^ A(06) = 0.01 GeV^ and A(08) = 0.01 GeV^ one obtains 
the relative errors 6Mo{9 GeV^) = 0.8% and 6Mi{9 GeV^) = 7.9%, where the latter is far 
beyond the accuracy needed to improve the pure data results on Aahad(^z) ^i-nd a^'^'^. 



Theoretical Uncertainties 

Looking at Eq. it is instructive to subdivide the discussion of theoretical uncertainties 
into three classes: 

{i) The perturbative prediction. The estimation of theoretical errors of the perturbative 
series is strongly linked to its truncation at finite order in a,. Due to the incomplete 
resummation of higher order terms, a non-vanishing dependence on the choice of 
the renormalization scheme (RS) and the renormalization scale is left. Furthermore, 
one has to worry whether the missing four-loop order contribution ds{as/7^Y gives 
rise to large corrections to the perturbative series. On the other hand, these are 
problems to which any measurement of the strong coupling constant is confronted 
with, while their impact decreases with increasing energy scale. The error on the 
input parameter a, therefore reflects the theoretical uncertainty of the perturbative 
expansion in powers of a,. A very robust a, measurement is obtained from the global 
electroweak fit performed at the Z-boson mass where uncertainties from perturbative 
QCD are reduced. The value found is as(M|) = 0.120±0.003 []19|. A second precise 



as measurement is obtained from the fit of the OPE to the hadronic width of the 
r lepton and to spectral moments , where the nonperturbative contribution was 
found to be lower than 1%. Additional tests in which the mass scale was reduced 
down to 1 GeV proved the excellent stability of the a, determination. The value 
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recently reported by the ALEPH Collaboration is a.(M|) = 0.1202 ± 0.0026 |2g . 
The consistency of the above values using quite different approaches at various 
mass scales is remarkable and supports QCD as the theory of strong interactions. 
Both measurements are almost uncorrelated so that we obtain the weighted average 
q;,(M|) = 0.1201 ± 0.0020, used in the following analysis. 

Although contained in the above uncertainty of a^, we add the total difference be- 
tween the results obtained using contour-improved fixed-order perturbation theory 



(FOPTci) and FOPT (see comments in Ref. |20|) as systematic error. 

Due to the truncation of the perturbative series, the arbitrary choice of the renor- 
malization scale /i leaves an ambiguity. When setting /x^ ^ s, additional logarithms 
enter the D-function (^. We evaluate the corresponding uncertainty by varying 
s<^^ < (3/2)s. 

{ii) The quark mass correction. Since a theoretical evaluation of the integrals and 
( P^ is only applied far from quark thresholds, quark mass corrections are small. 
We use the following settings: 

mu,d = , 
m,(l GeV) = (0.20 ± 0.07) GeV , 
m,(m,) = (1.3 ±0.2) GeV , 
mbimb) = (4.2 ± 0.2) GeV , 
mt{mt) = (176 ± 6) GeV . 

[in) The nonperturbative contribution. Using n = 1 only, the functions (|TT|) receive 
direct, i.e., non-suppressed contributions from the dimension D = 4 terms. Asso- 
ciated nonperturbative parameters are the vacuum expectation values of the gluon 
and the quark condensates. While the latter can be obtained from the PCAC rela- 
tion (|10D , for which a 20% uncertainty is assumed, the gluon condensate cannot be 
fixed theoretically. There fortunately exist experimental determinations using sim- 
ilar finite-energy sum rule techniques which are almost independent from the data 
used in this analysis: a fit using the r vector plus axial-vector hadronic width and 
spectral moments yields {{as/T^)GG) = (0.001 ± 0.015) GeV^ while a moment 
analysis using cc resonances results in {{as/'n')GG) = (0.017 ± 0.004) GeV^ EI 



An adjustment on e~^e data showed consistent results 0. Following the above 
estimates, we use a gluon condensate of 

{{as/n)GG) = (0.015 ± 0.020) GeV^ (15) 

in this analysis. 

Another tiny source of uncertainty is the error on the Z-boson mass. 



Low-Energy Results 

Due to the suppression of nonperturbative contributions in powers of the energy scale s, 
the critical domain where nonperturbative effects may give residual contributions to R{s) 
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is the low-energy regime with three active quark flavours. However, we have shown in 
Ref. [01 that at the scale ^/so = 1.8 GeV global duality holds and the small nonperturba- 
tive effects are described by the OPE. Up to this energy, R{s) is obtained from the sum 
of the hadronic cross sections exclusively measured in the occurring final states. 

The data analysis follows exactly the line of Ref. |1]], In addition to the e^e~ annihi- 
lation data we use spectral functions from r decays into two- and four final state pions 
measured by the ALEPH Collaboration |T^. As described in Ref. [Q, corrections to the 
charged width have to be applied to account for small CVC-violating effects. The 
magnitude of the width difference, (Tp± — rpo)/rp = (2.8 ± 3.9) x 10"'^ , translated into 
^had AQ;had(Af|) is evaluated using a parametrization of the p line shape based on 
vector resonances ITSII. One thus obtains the additive corrections 



fj, 



5a4;,(m|) 



-(1.3 ± 2.0) X 10"^° 
-(0.09 ±0.12) X 10"^ 



(16) 



for the T~ —>■ tt'tt^i't- spectral function which is applied in the present (and former) anal- 
ysis. Corrections for the higher mass resonances p(1450), p(1700) are expected to be 
negligible. Extensive studies have been performed in Ref. [0 in order to bound unmea- 
sured modes, such as KK with pions or the 7r+7r~47r° final states, via isospin constraints. 
We bring attention to the straightforward and statistically well-defined averaging proce- 
dure and error propagation used in this paper as in the preceeding ones, which takes into 
account full systematic correlations between the cross section measurements. All techni- 
cal details concerning the data analysis and the integration method used are presented in 
Ref. i. 

The experimental determination of the spectral moments in the first integral of Eq. (H) 
is performed as the sum over the respective moments of all exclusively measured e^e~ 
final states completed by r data. 

We fit the following equations {^/so = 1.8 GeV): 



so 



4m2. 



Piq[s 



a2(0) / ds 

|s|=so 



,had 



'^/i, [2m^, ^] 



Pil{s,) - Pi:>{s) D^Us) , (^ = 1, . . . , 9) (17) 



ds 



+ 



|s| = 



2\ Pl,o[^ 



s{s - Mz 
a(0)M| f ds r (,) 



Aahad(M|)[2m^, 



A\o(^5o) - Pi'lis) D^,,{s) , (i = 10, . . . , 18)(18) 



(19) 
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Mom. 


Cl 


Data 


A 


Theory 


A 


Sum 


A 


a 


(01) 





635.05 


7.36 


.00 


.00 


635.05 


7.36 


.00 


(02) 


.0005 


601.49 


6.93 


33.47 


.81 


634.96 


6.99 


.01 


(03) 


.001 


567.56 


6.51 


66.95 


1.62 


634.51 


6.73 


.07 


(04) 


.002 


500.18 


5.81 


133.89 


3.23 


634.07 


6.65 


.12 


(05) 


.003 


432.50 


5.27 


200.84 


4.85 


633.34 


7.17 


.19 


(06) 


.004 


364.98 


4.95 


267.78 


6.47 


632.77 


8.15 


.23 


(07) 


.005 


297.36 


4.90 


334.73 


8.08 


632.08 


9.46 


.27 


(08) 


.006 


229.78 


5.12 


401.67 


9.70 


631.45 


10.99 


.30 


(09) 


.007 


162.23 


5.59 


468.62 


11.32 


630.85 


12.64 


.31 


(10) 





56.77 


1.06 


.00 


.00 


56.77 


1.06 


.00 


(11) 


-.00005 


44.66 


.86 


11.98 


.29 


56.64 


.90 


.12 


(12) 


-.00007 


39.80 


.78 


16.78 


.41 


56.58 


.88 


.17 


(13) 


-.00009 


34.93 


.70 


21.57 


.52 


56.49 


.87 


.24 


(14) 


-.00011 


30.08 


.62 


26.36 


.64 


56.44 


.89 


.27 


(15) 


-.00013 


25.24 


.55 


31.15 


.75 


56.40 


.93 


.28 


(16) 


-.00015 


20.46 


.48 


35.95 


.87 


56.40 


.99 


.27 


(17) 


-.00017 


15.57 


.41 


40.74 


.99 


56.31 


1.07 


.31 


(18) 


-.00019 


10.68 


.35 


45.53 


1.10 


56.21 


1.16 


.35 



Table 1: Polynomial moments: The horizontal lines separate the moments into Eqs. ([T7|) 
and (0). The first column numbers the moments and the second defines the coefficients 
corresponding to Eq. (11). The following columns give the l.h.s integrals of Eqs. ( [TtD and 
(0), which are the experimental results, and the r.h.s integrals, i.e., the theoretical values 
after fitting. The sums of both experimental and theoretical integrals give the values for 



-had 



'^M2m^.i-8 Gev] ^ 10^° ^'^^ Aahad(M|)[2m^^ 1.8 GeV] X 10"^, Corresponding to the respective 
choice of polynomials. The last column shows the difference between the pure data results 
(given in the first and the 10th line, respectively, where the polynomial vanishes) and the 
sum in units of one standard deviation. 



where the polynomials Pi Q are defined in Eq. (|TT]) . Table |I| shows the measured polynomial 
moments (l.h.s. of Eqs. (|l^ and ([l^)), together with the fitted theoretical integrals (on 
the r.h.s. of Eqs. ( [TtI) and ([l8|)) and the sum of both integrals for all 18 polynomials. 
The correlation matrix including experimental and theoretical correlations of the moments 
defined in Table ^ is given in Table ^. These correlations have been inserted into the 
minimization used to fit the free varying parameters for which we obtain: 

^M2m^, 1.8 GeV] = (634.3 ± 5.6cxp ± 2.1thco) X 10 ^° , 
AQ;had(^z)[2m^, 1.8 GeV] = (56.53 ± 0.73exp ± 0.39thco) x 10"^ , (20) 

yielding a minimum = 1-0 for 16 degrees of freedom. The correlation between the 
fitted parameters (|20| ) amounts to 69%. 

The improvement in precision provided by this method can be easily understood from 
Table 0. As the chosen polynomial approaches better the actual kernel, the uncertainty 
from the integral on data is reduced, while the corresponding uncertainty from theory 
increases. Some optimum is reached for some particular choice of polynomial, providing 
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Mom. (01) (02) (03) (04) (05) (06) (07) (08) (09) (10) (11) (12) (13) (14) (15) (16) (17) (18) 



(1) 1 


.97 .94 


.82 


.65 


.46 


.31 


.19 


.10 


.75 


.71 


.66 


.59 


.51 


.42 


.34 


.25 


.18 


(2) - 


1 .97 


.88 


.73 


.57 


.43 


.31 


.23 


.71 


.71 


.68 


.63 


.57 


.50 


.42 


.35 


.29 


(3) - 


1 


.93 


.82 


.68 


.55 


.44 


.36 


.66 


.70 


.69 


.66 


.62 


.57 


.51 


.45 


.39 


(4) - 


- 


1 


.94 


.86 


.76 


.68 


.61 


.50 


.63 


.66 


.69 


.69 


.68 


.65 


.62 


.59 


(5) - 


- 


- 


1 


.95 


.90 


.85 


.80 


.31 


.51 


.59 


.65 


.70 


.73 


.74 


.73 


.73 


(6) - 








1 


.97 


.94 


.91 


.13 


.38 


.49 


.58 


.66 


.72 


.76 


.79 


.80 


(7) - 










1 


.97 


.96 


.00 


.28 


.40 


.52 


.62 


.70 


.76 


.80 


.82 


(8) - 












1 


.98 


-.10 


.19 


.33 


.45 


.57 


.67 


.74 


.79 


.83 


(9) - 


- 


- 


- 


- 


- 


- 


1 


-.17 


.13 


.27 


.40 


.53 


.63 


.72 


.78 


.82 


(10) - 
















1 


.93 


.86 


.78 


.67 


.56 


.44 


.34 


.24 


(11) - 


















1 


.97 


.93 


.87 


.79 


.70 


.62 


.54 


(12) - 




















1 


.97 


.93 


.87 


.80 


.73 


.66 


(13) - 
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.97 


.93 


.88 


.83 


.77 


(14) - 
























1 


.97 


.94 


.90 


.86 


(15) - 


























1 


.97 


.95 


.92 


(16) - 




























1 


.97 


.96 


(17) - 






























1 


.98 


(18) - 
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Table 2: Total experimental and theoretical correlations between the spectral moments used 
in the combined fit. 

the best constraint on the aj^^*^ and AQ;had(Af|) determinations. It should be noted that 
the experimental precision is not uniform throughout the energy range considered: for 
a/s ~ 0.8 GeV, R{s) is dominated by the p spectral function, well measured in both r 
decays and e~^e~ annihilation, while at larger s values the uncertainties from some poorer 
measurements such as e^e~ KKvrvr or 6tt take over. We clearly disagree with the 
conclusions reached in Ref. stating that the final uncertainty in the determination 
of Aahad(Af|) becomes insensitive to the experimental errors. The experimental uncer- 
tainty can hardly be reduced to nothing even if the integral on data becomes very small. 
A vanishing error is only achieved if the polynomial approximates well the kernel, but 
the theoretical uncertainty becomes prohibitive in this case. This trend is illustrated in 
Table |l|. 

As described before we may also use Eq. ([T3|) to constrain theoretically the low-energy 
dispersion relations (^) and (^). Here different moments are defined by varying values for 

and Ap. Unfortunately such moments are almost degenerate, i.e., they have very large 
correlations so that a combined fit as performed in the case of the polynomials does not 
make much sense. We therefore optimize the choice of and Ap in order to minimize 
the resulting errors on aj^|*2m,, ^ and Aahad(M|)[2rr»^, ^ setting again ^ = 1.8 GeV. 
The corresponding results are given in Table 0. An analysis of the theoretical uncertain- 
ties shows that at \/Q^ = 3 GeV nonperturbative contributions to the D-function are 
negligible. The theoretical errors given in Table |^ originate mainly from uncertainties of 
the massless perturbative series in Eq. (^), i.e., the error on as and the variation of fi. The 
weighted averages of pure data and the theory- improved results given take into account 
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the correlations between both values. These values agree with those from the combined 
fit of the polynomial moments. The most precise of the two evaluations are retained for 
the final results. 



Results on the cc Threshold 

The cc threshold region involves the narrow resonances, J/ipllS), ip{2S) and ijj{3770), 
parametrized using relativistic Breit-Wigner formulae |2^, 0, as well as the broad reso- 
nance and continuum region starting with the opening of the e"*"e~ —>■ DD mode at about 
3.7 GeV. A compilation of the data used can be found in Ref. [0]. In this energy region we 
do not use the polynomial approach for which the precision is limited by the uncertainty 
on the c-quark mass and the threshold behaviour of the correlator. The first problem 
is less severe with the dispersion relation approach while the second one is inexistent. 
Equation (pTsD now reads: 



S2 



ds R{s) 



fis) 



+Af Id,{Q^)-Q^ J ds 



+ E 

V'=15,25,3770 



ds R^{s) 



fis) 



{s + Q 



2^2 



S2 

+ J ds 



QCD/ 

uds \ 



si 



{s + Q 



2\2 



AfQ' 



(21) 



where we choose for the continuum the integration ranges ^/si = 3.7 GeV and ^/s2 = 
5 GeV, and set F = a^'^f^^i^^ ^ and F = Aahad(M|)[^][^, while /(s) denotes the 
corresponding integration kernels, respectively. Table ^ gives the results of the individual 
terms in Eq. (|2T|) again after optimizing the choice of and the normalization Ap. The 
quoted error on Dc{Q'^) is dominated by the uncertainty on the c-quark mass. 



Results for Aahad(^z) '^Ji^^ 

In the previous sections we reevaluated the dispersion integrals yielding the hadronic 
contributions to and to a{M^) for the low-energy and cc threshold regions. According 
to Ref. H] (see also Ref. [P3); employ perturbative QCD using the formulae of Ref. P 



for the cross section ratio R{s) in other regions. For the crossing of the bb threshold, 
we assume nonperturbative effects to be negligible and use perturbative QCD in the 



context of global duality (see, e.g., Ref. |^) to evaluate the hadronic contributions. The 
theoretical error is then dominated by the uncertainty on the 6-quark production threshold 
2Mh, with Mjy being the pole mass of the 6-quark. 

Table |^ shows the experimental and theoretical evaluations of Aa^adiM^), a^^'^ and 
^had ^YiQ respective energy regimes^. Experimental errors between different lines are 
assumed to be uncorrelated, whereas theoretical errors but those from cc and bb thresholds 
which are quark mass dominated are added linearily. 



^ The evaluation of a^'^'^ follows the same procedure as a^^'^. 
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Aahad(iWz)[2m,, 1.8 GeV] X 10 


"m, [2m^. 1-8 GeV] 




3 GeV 


3 GeV 


Af 


-0.0003 


0.009 


P^P (2m^ - 1.8 GeV) 


44.67 ±0.70 


534.1 ±6.2 


42.18 ±0.23 


353.5 ±2.0 


(1-8 GeV-oo) 


-30.68 ±0.21 


-257.1 ± 1.8 


Total 


56.17 ±0.70exp±0.22theo 


630.5 ±6.2exp±1.9theo 


Data only 


56.8 ± 1.1(1) 


635.1 ± 7.4(1) 


Correlation 


95% 


95% 


Average 


56.36 ±0.70exp±0.18thco 


632.5 ± 6.2cxp ± 1.6thoo 






Aahad(M|)[^][3.7, 5 GeV] X 10^ 


had y 1 nlO 
"m, [V][3.7, 5 GcV] X iU 




15 GeV 


15 GeV 


Af 


-0.0009 


0.0015 


r^'P (3.7-5 GeV) 


7.02 ±0.45 


2.94 ±0.19 


J^^P (V^(l^,2^,3770)) 


6.25 ±0.37 


6.12 ±0.35 


79.95 ± 0.83 


37.22 ± 0.38 


-Jf ^° (5 GeV - oo) 


-74.01 ±0.31 


-34.55 ±0.15 


Jfr (3.7 - 5 GeV) 


5.32 ±0.03 


2.47 ±0.02 


Total 


24.53 ± 0.58exp ± 0.88theo 


14.20 ± 0.40exp ± 0.40theo 


Data only (3.7-5 GeV) 


15.80 ± 1.00(2) 


6.93 ±0.44(2) 


Data only ^(15,2^,3770) 


9.24 ±0.68(3) 


7.51 ± 0.44(3) 


Correlation 


55% 


71% 


Average 


24.75 ± 0.84exp ± 0.50thco 


14.31 ± 0.50exp ± 0.21theo 



^ The slight modification compared to om- previous values of AQ;had(Ai^|)[2m„, i.s GcV] = 
56.9 X 10"'' and 

^^^fzm^ 1 8 GcV] ~ 636.5 X 10 1" are due to a reevaluation of the contri- 
butions from T data 

^ A reevaluation of the treatment of the experimental errors on the contribution from charm 
threshold led to a significant reduction of the very conservative uncertainties given in our 
previous paper Q 

The changes in the resonance contributions compared to Ref. is due to the multiplicative 
correction (q!/q;(M^))^, erroneously not applied before [ p3[ 

Table 3: The optimized choices of the parameters and Af, the solutions of the corre- 
sponding integrals ('T") in Eq. (|I3|) (upper table) and Eq. (^) (lower table) as well as the 
corresponding values for the Adler D-functions. Additionally given are the correlations 
between the theory-improved and the pure data results used to calculate the average of both. 
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Energy (GeV) Aahad(M|) x 10^ aj^^^ x lO^o x 10^^ 

(2m^ ^ l.S)uds 56.36 ± 0.70exp ± O.lStheo 634.3 ± 5.6exp ± 2.1theo 173.67 ± 1.7exp ± 0.6theo 

(1.8 - 3.700),rf, 24.53 ± 0.28thco 33.87 ± 0.46theo 8.13 ± O.lltheo 

V^(15 25 3770)e 24.75 ± 0.84exp ± 0.50theo 14.31 ± 0.50exp ± 0.21theo 3.41 ± 0.12exp ± 0.05theo 

(5-9.3)«rf,c 34.95 ± 0.29theo 6.87 ± O.lltheo 1.62 ± O.OStheo 

{92 - I2)udscb 15.70 ± 0.28theo 1.21 ± 0.05theo 0.28 ± 0.02theo 

[12 -^)udscb 120.68 ± 0.25theo 1.80 ± O.Oltheo 0.42 ± O.Oltheo 

(2mt - oo)t -0.69 ± 0.06thco ~ « 

(2m7r - oo)udscbt 276.3 ± l.loxp ± 1-lthco 692.4 ± 5.6oxp ± 2.6theo 187.5 ± 1.7cxp ± 0.7thoo 



Table 4: Contributions to Aahad(^z)j '^Ji^*^ ^'^^ '^e'^'^ /rom i/ie different energy regions. 
The subscripts in the first column give the quark flavours involved in the calculation. 



According to Table ^ the combination of the theoretical and experimental evaluations 
of the integrals (|^) and @ yields the final results 



Afthad (M|) = 


(276.3 ±l.lexp±l.ltheo) X 10-4 


a-\Mi) = 


128.933 ± O.OlSexp ± 0.015theo 


„had 


(692.4 ± 5.6exp ± 2.6theo) x IQ-^o 


«r = 


(11 659 159.6 ± 5.6cxp ± 3.7theo) x 10-^° 



and a^^'^ = (187.5 ± 1.7exp ± 0.7theo) x IQ-^^ for the leading order hadronic contribu- 
tion to Og. The total value includes an additional contribution from non- leading 
order hadronic vacuum polarization summarized in Refs. p6|, |l]| to be a^^'^[{a/n)^] = 
(—10.0 ± 0.6) X 10-^°. Also the light-by-light scattering (LBLS) contribution has re- 
cently been reevaluated to be a|j^'^[LBLS] = (-7.9 ± 1.5) x 10"^° Together with 
the value aJ^^[LBLS] = (-9.2 ± 3.2) x 10-^° ||2|, we use the average (aJ^'^[LBLS]) = 
(—8.5 ± 2.5) X lO"^*^ so that the total higher order hadronic correction amounts to 
a^^'^^lia/nf + LBLS] = (-18.5 ± 2.6) x 10"^°. Figures | and | show a compilation of 
published results for the hadronic contributions to a(M|) and a^. Some authors give the 
hadronic contribution for the five light quarks only and add the top quark part separately. 
This has been corrected for in Fig. |^. 
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Figure 2: Comparison o/Aahad(^z) evaluations. The values are taken from Refs. 
2|, BOl pll H, [I, I, m, 1^ and from this work. 
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Figure 3: Comparison of a^'^ evaluations. The values are taken from Refs. [Q, 
22| , ^ |l], H and from this work. 
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Conclusions 



We have reevaluated the hadronic vacuum polarization contribution to the running of the 
QED fine structure constant, a{s), at s = M| and to the anomalous magnetic moment 
of the muon, a^. We employed perturbative and nonperturbative QCD in the frame- 
work of the Operator Product Expansion in order to extend the energy regime where 
theoretical predictions are reliable. Based on analyticity, we constrained the low-energy 
and cc threshold region theoretically using finite energy sum rules and dispersion rela- 
tions. The standard evaluation using data from e~^e~ annihilation and r decays at low 
energies and near quark production thresholds is therefore improved. Our results are 
Aahad(M|) = (276.3 ± 1.6) x 10-^ propagating ^-^(O) to a-i(M|) = (128.933 ± 0.021), 
and al""^ = (692.4 ± 6.2) x 10^° which yields the Standard Model prediction af^ = 
(11659159.5 ± 6.7) x 10-^°. For the electron we found a^''^ = (187.5 ± 1.8) x 10"^^ 
These results have direct implications for on-going experimental programs. On one hand, 
the precision on a(M|) is now such that it does not limit anymore the adjustment of 
the Higgs mass from accurate experimental determinations of sin^6'w 0. On the other 
hand, the gain in accuracy for a^^'^ is even more rewarding as it will permit to exploit the 
foreseen precision increase of the measurement [^] in order to achieve a significant 
determination of the contribution a^'^^^. 
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